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摘要 :【 目 的] 为 了 系统 了 解 中 华 蜜蜂 Apis cerana cerana 蜂王 和 雄 蜂 转录 组 特征 ,丰富 蜜蜂 转录 组 数据 信息 。[ 方 
法 ] 本 人 研究 利用 高 通 量 测序 的 方法 分 别 检测 中 华 蜜蜂 蜂王 和 雄 蜂 刚 出 房 .性 成 熟 时 期 以 及 性 成 熟 期 蜂王 生殖 系统 
和 雄 蜂 生 殖 系 统 之 间 转 录 组 表达 差异 。[【 结果 】 经 过 测序 获得 质量 值 不 低 于 20 的 碱 基 比 例 (Q20) 均 高 于 90% ;所 
有 reads 组 装 成 90 839 个 unigenes ,平均 长 度 1 549 bp; 基 于 5 个 数据 库 ( NR ,Swiss-Prot,GO,COG 和 KEGG ) 进行 比 
对 ,共有 45 112 个 unigenes 被 注释 。 差 异 基因 表达 分 析 发 现 , 与 刚 出 房 时 相 比 ,性 成 熟 的 蜂王 和 雄 蜂 均 在 表皮 和 蛋 
白 .细胞 色素 P450 .气味 结合 蛋白 等 家 族 基因 上 存在 显著 差异 表达 ,而 且 这 些 差异 表达 基因 与 蜜蜂 生长 发 育 和 性 成 
熟 过 程 中 蜜蜂 骨骼 发 育 .生殖 系统 发 育 嗅觉 发 育 等 方面 有 关 ; 性 成 熟 蜂王 与 性 成 熟 雄 峰之 间 以 及 它们 生殖 系统 之 
间 在 气味 结合 蛋白 基因 方面 存在 显著 差异 。[【 结论 ] 结 果 表 明 , 中 华 蜜蜂 在 性 成 熟 过 程 中 ,体内 大 量 基因 的 表达 发 
生 了 变化 。 这 些 结果 揭示 了 中 华 蜜蜂 性 成 熟 发 育 的 整体 基因 表达 特征 ,在 得 到 大 量 转录 组 unigene 序列 的 同时 , 获 
得 了 一 批 与 蜜蜂 性 成 熟 有 关 的 基因 序列 ,为 深入 开展 中 华 蜜蜂 生长 发 育 与 繁殖 研究 提供 了 丰富 的 数据 资源 。 
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Abstract: [Aim) This study aims to understand the transcriptome characteristics of queens and drones of 
the Chinese honeybee, Apis cerana cerana, and enrich the transcriptome data. [ Methods 】 The 
transcriptomes of queens and drones newly emerged and sexually matured were constructed and analyzed 
by RNA-seq. [Results) The results showed that the proportion of nucleotides with quality value larger 
than 20 in reads ( Q20) were more than 90% and all the reads were assembled into 90 839 uniques with 
a mean length of 1 549 bp. Based on sequence similarity search against five public databases (NR, 
Swiss-Prot, GO, COG and KEGG), a total number of 45 112 unigenes were annotated. Using these 
transcriptome data as the reference sequences, we found that the genes encoding cuticular proteins/ 
apidermins, CYP and odorant binding proteins were differentially expressed between the sexually matured 
bees and the newly emerged bees which are related to the development of body, reproductive system and 
olfaction system. Meanwhile, the odorant binding protein genes were significantly differentially expressed 
between the sexually matured queens and drones as well as their reproductive system. 【 Conclusion 】The 
results indicated that the expression levels of a large number of genes changed during maturing of A. cerana 
cerana bees. This study made an insight into the characteristics of the development in adult queen and 
drone, and a lot of transcript sequences with important function were acquired for future gene expression or 
regulation research relative to the growth, development and reproduction in A. cerana cerana. 
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中 华 蜜 蜂 Apis cerana cerana 是 我 国 特 有 的 宝贵 Apis mellifera 相 - 中 华 蜜蜂 具有 嗅觉 灵敏 ,善于 利 
蜜蜂 资源 ,也 是 东方 蜜蜂 指名 亚 种 。 与 西方 蜜蜂 ”用 零星 蜜 粉 源 ,飞行 临界 温度 低 , 有 利于 冬季 部 分 植 
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物 授 粉 和 冬 密 生产 , 蜂 群 的 抗 螨 性 能 强 等 优势 ( 曾 
志 将 ,2009) 。 然 而 ,由 于 杀 虫 剂 和 除草 剂 等 在 现代 
农业 中 的 广泛 应 用 ,给 蜜蜂 健康 养殖 带 来 了 巨大 威 
胁 (Frazier et al.，2008; 代 平 礼 等 ,2009; Henry et 
al., 2012; Whitehorn et al., 2012)。 男 外 ,由 于 西方 
蜜蜂 的 进入 ,导致 中 华 蜜蜂 种 群 数 量 与 分 布 区 域 明 
显 缩小 ( 季 荣 等 ,2003; 杨 冠 焊 ,2005 ) ,因而 ,合理 利 
用 和 保护 我 国 这 特有 的 蜜蜂 资源 势 在 必 行 。 

蜜蜂 为 了 提高 蜂 群 的 延续 性 和 适应 性 ,已 进化 
出 其 特殊 的 竞争 交尾 机 制 , 即 婚 飞 行为 。 蜂 王 和 雄 
蜂 只 有 到 达 性 成 熟 , 才 会 选择 合适 的 机 会 进行 婚 飞 ， 
这 也 说 明 在 性 成 熟 发 育 过 程 中 ,蜂王 和 雄 蜂 体内 发 
生 了 明显 的 生理 变化 ,而 且 蜜蜂 在 不 同 发 育 阶 段 ,其 
体内 相关 基因 表达 调控 也 存在 明显 区 别 ,并 已 经 在 
蜜蜂 生理 生化 和 分 子 水 平 上 得 到 证 实 ( Kocher et 
al., 2008,2010; Colonello-Frattini and Hartfelder, 
2009; Behura and Whitfield, 2010; Chen et al., 
2010;Fang et cl.，2012; 吴 小 波 等 ,2012; 曾 晶 等 ， 
2012 ; 石 元 元 等 ,2012; 刘 芳 ,2012) 。 性 成 熟 处 女 蜂 
王 在 进行 婚 飞 过 程 中 ,其 体内 发 生 了 大 量 的 基因 表 
达 差 异 ( 吴 小 波 等 ,2013a， es et al.，2014 ) 。 
研究 还 发 现 ,西方 塞 蜂 会 干扰 中 华 蜜蜂 处 女 蜂王 交 
配 ， bo i 分 泌 的 性 引诱 

息 素 不 但 会 吸引 中 华 蜜蜂 的 雄 蜂 来 交尾 ,而且 还 
会 吸引 大 量 西方 蜜蜂 的 雄 蜂 ,从 而 使 中 华 蜜蜂 蜂王 
交尾 成 功率 大 幅度 下 降 ( 王 启发 等 ,2003 ) 。 蜂 王 和 
雄 蜂 作为 两 种 性 别 的 个 体 ,其 性 成 熟 时 期 体内 的 生 
理 状况 也 不 一 样 ,但 对 其 差异 及 相关 调控 机 理 尚未 
深入 研究 。 

近年 来 , 随 着 高 通 量 测序 技术 的 不 断 发 展 与 
完善 ,昆虫 的 转录 组 学 得 到 快速 发 展 ( 张 棋 鹿 等 ， 
2013) ,而且 蜜 蜂 转 录 组 研究 取得 了 阶段 性 研究 
成 果 , 但 关于 中 华 蜜蜂 转录 组 研究 报道 较 少 , 仪 对 中 
华 蜜蜂 整体 转录 情况 进行 了 初步 探索 ( Wang et al.， 
2012) ,缺少 对 中 华 蜜蜂 不 同 发 育 阶段 的 转录 组 
表达 信息 的 了 解 , 尤 其 是 性 成 熟 发 育 阶 段 。 本 研究 
拟 通过 转录 组 测序 , 全面 探讨 中 华 蜜蜂 蜂王 和 雄 
蜂 在 性 成 熟 时 与 刚 出 房 时 的 转录 组 基本 信息 ,并 初 
步 分 析 了 性 成 熟 时 期 蜂王 生殖 系统 与 雄 蜂 生殖 系统 
的 转录 组 信息 ,了 解 中 华 蜜蜂 在 性 成 熟 过 程 中 特定 
基因 的 表达 情况 ,为 深入 研究 中 华 蜜蜂 性 成 熟 的 发 
育 机 制 、 主 要 相关 基因 克隆 以 及 功能 分 析 等 黄 定 
基础 。 












































































































































1 材料 与 方法 


1.1 供 试 昆虫 

实验 蜂 群 为 江西 农业 大 学 蜜蜂 研究 所 饲养 的 中 
华 蜜蜂 ,其 蜂 群 各 含有 1 头 同 批 次 且 交 尾 成 功 的 新 
蜂王 。 冤 蜂 自由 出 党 采集 ,不 进行 人 工 饲 喂 。 
1.2 样品 的 采集 

蜂王 :按照 标准 的 人 工 育 王 方 法 培育 中 华 蜜蜂 
蜂王 ( 曾 志 将 ,2009) , 待 中 华 蜜蜂 蜂王 出 台 后 ,立即 
采集 刚 出 房 的 蜂王 ， 迅速 用 解剖 剪 从 蜂王 侧面 解剖 
并 剔除 蜂王 体内 的 消化 道 , 放 和 人 液 氮 中 保存 , 待 提 
取 RNA。 

按照 同样 的 方法 培育 中 华 蜜蜂 蜂王 ,将 新 出 房 
Pe ee 处 女 蜂 王 出 房 
12 d 后 ,采集 性 成 熟 处女 蜂 王 ,并 剔除 蜂王 体内 的 
消化 道 , 放 入 液 氮 中 保存 , 待 提 取 RNA。 
雄 蜂 ; 用 隔 王 板 将 已 开始 产 雄 蜂 卵 的 蜂王 控制 
在 带 有 雄 蜂 站 房 的 巢 脾 中 产 卵 ,24 h 之 后 ,移出 蜂 
王 ,并 将 该 梨 脾 放 在 中 间 继 续 孵 化 发 育 。 待 雄 蜂 即 
将 出 房 时 , 拌 落 集 脾 中 的 工蜂 ,将 该 巢 脾 放 入 恒温 恒 
温 箱 (温度 35%C ,RH 95% ,无 光照 ) 中 ,让 雄 蜂 羽化 
出 房 ,并 采集 刚 出 房 的 雄 蜂 。 同 时 ,对 部 分 刚 出 房 雄 
蜂 进 行 标记 , 放 入 原 实 验 种 群 进一步 发 育 ,出 房 16 d 
后 ,采集 这 些 已 标记 的 雄 蜂 ,用 解剖 剪 从 雄 蜂 侧面 解 
齐 并 剔除 雄 蜂 体 内 的 消化 道 , 放 入 液 氮 中 保存 , 待 提 
取 RNA。 

蜜蜂 生殖 系统 :通过 控制 蜂王 产 雄 蜂 卯 和 人 工 
培育 蜂王 ,获取 性 成 熟 的 蜂王 和 雄 蜂 ,迅速 用 解剖 剪 
从 蜜蜂 侧面 解剖 并 剔除 体内 的 消化 道 ,取出 生殖 系 
统 , 采 集 蜂 王 生 殖 系 统 和 雄 蜂 生殖 系统 样品 , 放 入 液 
毛 中 保存 , 待 提取 RNA。 
1.3 样品 RNA 提取 

从 每 组 样品 中 随机 挑选 4 个 样品 ,按照 标准 的 
方法 分 别提 取 蜂 王 和 雄 蜂 刚 出 房 时 期 .性 成 熟 时 期 
以 及 生殖 系统 样品 的 总 RNA, 并 进行 RNA 质量 检 
测 。 待 样品 检验 合格 后 ,将 每 组 试验 组 的 两 个 样品 
等 量 混合 RNA ,各 制备 2 个 平行 样 。 
1.4 创建 cDNA 文库 和 转录 组 测序 

分 别 将 各 试验 样品 中 RNA 进行 cDNA 反 转 录 ， 
并 进行 纯化 等 处 理 , 然后 用 琢 脂 糖 凝 胶 电 泳 对 片段 
大 小 进行 选择 , 最 后 进行 PCR 扩 增 ,获得 最 终 的 
cDNA 文库 。 通 过 Solexa RNA 的 paired-end 测序 方 
法 对 各 试验 样品 进行 测序 ,并 对 测序 获得 原始 数据 
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经 过 去 除 低 质量 片段 ,切除 reads 中 含 N 部 分 序列 
后 ,获取 干净 序列 。 具 体 参 照 Wang 等 (2012 ) 研究 
方法 。 
1.5 中 华 蜜蜂 性 成 熟 期 差异 表达 基因 分 析 

将 试验 获取 的 Clean reads 与 Wang 等 (2012 ) 测 
序 获得 的 中 华 蜜 蜂 转 录 组 数据 进行 混 拼 ,然后 以 此 
为 参照 ,利用 RPKM 法 计算 每 个 unigene 的 表达 量 
(Mortazavi et al.,，2008 ) 。 在 此 基础 上 ,参照 Audic 
和 Claverie(1997 ) 等 描述 的 方法 ,对 蜂王 、 雄 蜂 性 成 
熟 过 程 中 以 及 性 成 熟 期 生殖 系统 差异 表达 基因 进行 
筛选 ,并 对 这 些 差异 基因 在 CO( Gene Ontology ) 数据 
库 进 行 基因 功能 富 集 分 类 , 在 KEGG ( Kyoto 
Encyclopedia of Genes and Genomes ) 数据 库 进 行 代 
谢 途 径 富 集 分 析 。 
1.6 ”qRT-PCR 验证 

为 了 验证 试验 结果 的 可 靠 性 ,分 别 重新 采集 上 
述 样品 ,经 过 剔除 消化 道 后 ,按照 1.2 和 1.3 节 的 方 
法 制备 qRT-PCR 验证 试验 样品 。 从 比 对 结果 中 选 
择 了 5 个 共同 差异 基因 ,在 性 成 熟 蜂 王 与 刚 出 房 蜂 
王 以 及 性 成 熟 蜂 王 生殖 系统 与 性 成 熟 雄 蜂 生 殖 系 统 
样品 中 进行 试验 验证 ,并 进行 统计 分 析 。 引 物 设计 
如 表 1。 


















































2 结果 与 分 析 


2.1 中 华 蜜蜂 转录 组 Jumina 测序 和 序列 组 装 

对 中 华 蜜蜂 羽化 和 性 成 熟 时 期 的 蜂王 和 雄 蜂 样 
品 进 行 Humina 测序 ,将 测序 获得 的 原始 图 像 数 据 
经 basecalling 转 化 为 序列 数据 ,并 进行 过 滤 去 杂 











表 1 实时 荧光 定量 PCR 引物 
Table 1 qRT-PCR primers 
基因 引物 序列 (5 -3") 


Cene Primer sequence 


F: CCAACCATCGAAGAAGTAGAC 








HSP70Ab 
R: ATGACTGTGCAAAGTTTGCAG 
ed F: ATTGAAGGAGTTTGATGGCA 
R: AGAGATTCTCGAATTTGGCT 
pr F: CGATTGGTTGATTCACCTT 
R: GCAGCCATATATCCCATAGTAG 
BE F: ATTCCATTTCGTGGTCTCT 
R: TGAACTAAACGCAAAGTCTCA 
ee F: CTAGTTTATCATCGGATGGTGT 
SR R: CCTTCAATGCCGGTTTAC 
F: TTCATGCTGTTACTGCTACAC 
GAPDH 


R: GAAGGCCATACCAGTCAAT 





(包括 含有 带 接头 重复 的 和 测序 质量 很 低 的 
reads) ,得 到 干净 序列 。 从 表 2 中 可 以 看 出 ,经 过 测 
序 获取 到 的 干净 序列 数 在 43 220 752 ~ 51 791 872 
之 间 , 累计 长 度 介 于 4 322 075 200 ~5 179 187 200 
bp 之 间 。 测 序 质量 值 大 于 等 于 20( Q20 ) 的 碱 基 比 
例 均 高 于 94% ,GC 含量 36% ~41% 之 间 ( 表 2)。 
对 测序 数据 进行 组 装 后 分 析 ( 表 3) ,本 研究 获取 的 
每 个 样品 的 clean reads 能 组 装 成 contigs 数 介 于 88 
361 ~ 107 490 ,平均 长 度 396 ~432 bp 之 间 。Contigs 
的 N50 介 于 924 ~ 1 114 bp。 将 这 些 contigs 和 Wang 
等 (2012) 实验 原始 数据 混 拼 ,通过 双 末 端 链接 和 间 
院 填 充 最 终 组 装 成 90 839 个 unigenes ,平均 长 度 为 
1 549 bp ,unigenes 的 N50 是 3 156 bp。 对 最 后 组 装 
出 来 的 unigene 进行 长 度 分 布 统计 ,共有 41 536 个 
unigenes 长 度 超过 了 1 000 bp( 图 1)。 









































表 2 中 华 蜜蜂 转录 组 数据 分 析 


Table 2 Analysis of the transcriptome data of Apis cerana cerana 























样品 原始 序列 总 数 。 十 兆 序列 总 数 。 王 淮 核 帮 总 数 (n) 020 NN 百分比 6C 百 分 此 
Samples Total raw reads Total clean reads Total clean nucleotides N percentage GC percentage 
1 由 由 条 
刚 出 房 雄 蜂 58 304 712 48 814 742 4 881 474 200 96. 83 0.01 37599 
Newly emerged drone 
| TH s 
刚 出 房 蜂 王 59 197 522 51 403 028 5 140 302 800 97.38 0.01 37.16 
Newly emerged queen 
性 成 就 雄 峰 
性 成 熟 雄 蜂 58 470 882 51 407 440 5 140 744 000 97.58 0.01 38.35 
Sexually matured drone 
E 蜂 和 牛 殖 系统 
雄 峰 后 殖 系统 59 930 148 51 791 872 5 179 187 200 97.30 0.01 40.16 
Drone reproductive system 
王后 统 
蜂王 生殖 系统 61 966 372 45 534 872 4 553 487 200 94.96 0.04 36.16 
Queen reproductive system 
性 成 训 峰 干 
性 成 熟 蜂 王 58 121 396 43 220 752 4 322 075 200 94.98 0.04 37.02 
Sexually matured queen 
参考 组 
51 581 510 51 581 510 4 642 335 900 97.55 0.03 41.92 
Reference 














Q20: 质量 值 不 低 于 20 的 碱 基 比例 Proportion of nucleotides with the quality value larger than 20 in reads. 参考 组 为 Wang 等 (2012 ) 研究 中 华 蜜蜂 
转录 组 序列 结果 ;下 同 。The reference is the transcriptome of the Asian honey bee by Wang et al. (2012). The same below. 
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表 3 测序 序列 组 装 质量 
Table 3 Assembly quality of sequencing sequence 








样品 总 数 总 长 度 (nt) 平均 长 度 (nt) N50 
Samples Total number Total length Mean length (nt) 
| 
Contig 刚 出 房 雄 蜂 92 081 39 783 756 432 1 114 
Newly emerged drone 
1 
刚 出 房 蜂王 88 361 37 792 699 428 1 085 
Newly emerged queen 
性 成 就 雄 
性 成 熟 雁 蜂 105 089 42 245 210 402 999 
Sexually matured drone 
雁 峰 生殖 系 4 96 339 40 236 929 418 1 016 
Drone reproductive system 
蜂王 生殖 系统 107 490 42 611 757 396 924 
Queen reproductive system 
性 成 熟 
性 成 熟 蜂王 97 693 40 344 556 413 1 015 
Sexually matured queen 
参考 
参考 组 76 087 25 143 199 330 528 
Reference 
| i 
Unigene 刚 出 房 雄 蜂 52535 75 438 695 1 436 2 778 
Newly emerged drone 
1 
刚 出 房 蜂王 48 213 64 088 959 1 329 2 558 
Newly emerged queen 
性 成 就 雄 
性 成 熟 雁 蜂 55 689 77 987 467 1 400 2 716 
Sexually matured drone 
雄 峰 生殖 系统 54 413 74 298 541 1 365 2 587 
Drone reproductive system 
蜂王 生殖 系统 56 182 67 709 383 1 205 2 337 
Queen reproductive system 
性 成 熟 
性 成 熟 蜂王 52 180 66 632 174 1 277 2 459 
Sexually matured queen 
参考 
参考 组 36 972 27 233 479 737 1 124 
Reference 
合计 
90 839 140 715 277 1 549 3 156 
Total 


“把 组 装 出 的 contigs 或 scaffolds 从 大 到 小 排列 , 当 其 累计 长 度 刚 刚 超过 全 部 组 装 序列 总 长 度 50% 时 ,最 后 一 个 contig 或 scaffold 的 大 小 即 为 
N50 的 大 小 。The assembled contigs or scaffolds are arranged from large to small. When the cumulative length is over 50% of the total length of all 


assembly sequences, the size of the last contig or scaffold is N50. 
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图 1 Contigs 和 unigenes 从 头 组 装 的 长 度 分 布 


Fig.1 Length distribution of the de novo assembly for contigs and unigenes 
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2.2 Unigene 功能 注释 

将 unigene 分 别 注释 到 NR ,Swiss-Prot ,KECGC ， 
COG 和 GO 数据 库 中 ,并 分 别 对 注释 到 每 个 库 以 及 
所 有 注释 上 的 unigene 数目 进行 统计 。 总 共 获 得 
45 112 个 unigenes ( 占 49. 66% ) 具有 比 对 结果 ( 表 
4) 。 其 中 通过 NR ,Swiss-Prot,GO ,COG 和 KEGG 数 
据 库 比 对 到 的 结果 分 别 是 44 774(49.26% ) ,36 294 
(39. 95% ) ,20 850 (22.95% ) ,19 169 (21. 10%) 
和 33 448 (36. 82% ) 。 

表 4 Unigenes 的 注释 结果 
Table 4 Annotation of unigenes 


注释 unigenes 的 百分比 





数据 库 注释 unigenes 的 数目 





Percentage of 





Database Number of annotated unigenes ee 
NR 44 744 49.26 
Swissprot 36 294 39. 95 
COG 19 169 21.1 
GO 20 850 22.95 
KEGG 33 448 36. 82 
Total 45 112 49.66 


被 NR 注释 到 的 unigenes 基因 中 ,相似 性 高 于 
95% 基因 占 51.5% ,相似 性 介 于 80% ~95% 之 间 的 
基因 数 占 28.3% (图 2) 。 将 这 些 unigenes 比 对 到 各 
个 物种 中 ,分别 有 53.7% ,26.4% 和 8.4% 的 基因 比 
对 到 西方 蜜蜂 、 小 蜜蜂 中 和 熊 蜂 中 ,另外 还 有 7.6% 
的 未 知 基 因 ( 图 3)。 将 这 些 unigenes 在 COG 数据 
库 中 进行 分 类 , 共 分 25 类 (图 4) ,其 中 最 大 群体 是 
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2.3 ”中华 蜜蜂 不 同时 期 差异 表达 基因 分 析 


为 了 检测 中 华 蜜蜂 不 同时 期 基因 表达 差异 ,我 


Rl 


R 类 (一 般 功 能 预测 ) ,有 8 201 个 基因 ,其 次 是 上 L,] 
和 KK 类 (“复制 ,重组 和 修复 ”“ 翻译 ,核糖 体 结 构 和 
生物 起 源 ”" “转录 ”) ,分 别 对 应 3 748,3 340 和 
3 277 个 基因 。 进 一 步 将 注释 到 GO 数据 库 中 的 
20 850 个 unigenes 按照 生物 过 程 、 细 胞 组 成 分 子 功 
能 三 大 类 别 进行 分 类 ,其 中 按照 生物 过 程 分 成 23 个 
类 别 ,按照 细胞 组 成 分 成 18 个 类 别 ,按照 分 子 功能 
分 成 18 个 类 别 。 
相似 性 (9%) 
Similarity 
国 17%-40% 
园 40%-60% 
口 60%-80% 


口 80%-95% 
国 95 和 -100% 


7.30% 
0% 







.80% 


51.50% 


28.30% 


图 2 注释 到 的 unigenes 基因 相似 性 分 布 


Fig.2 Similarity distribution of the annotated unigenes 








O80% -7.60% 
pe 2 9 西方 蜜蜂 Apis mellifera 
2 a 小 蜜蜂 Apis florca 


o 能 蜂 Bombus bee 

o 首 若 切 叶 蜂 Megachile rotundata 
日 切 叶 凡 Acromyrmex echinatior 

日 印度 跳 蚊 Harpegnathos saltator 


日 其 他 0thers 
“53.70% 
图 3 注释 到 的 unigenes 基因 种 类 分 布 


Fig.3 Species distribution of the annotated unigenes 
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图 4 ” ”Unigenes 的 COG 分 类 
Fig.4 COG classification of the unigenes 


们 分 析 了 中 华 蜜蜂 蜂王 和 雄 蜂 在 性 成 熟 期 与 刚 出 房 
时 差异 表达 基因 ,并 对 同一 时 期 两 种 性 别 的 中 华 蜜 
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蜂 及 生殖 系统 进行 差异 基因 分 析 , 具体 差异 数量 如 
图 5 所 示 , 从 图 中 可 以 看 出 :蜂王 和 雄 蜂 在 性 成 熟 其 
上 调 表 达 的 基因 均 比 刚 出 房 时 要 多 ; 刚 出 房 的 雄 蜂 
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上 调 表达 的 基因 比 蜂王 多 ,但 性 成 熟 期 的 雄 蜂 上 调 
表达 的 基因 上 比 蜂 王 少 ;性 成 熟 期 蜂王 卵 梨 上 调 表 达 
的 基因 比 雄 蜂 精 集 要 多 。 








国 上 调 Up 
国 下 调 Down 





SD-ED 


EQ-ED ED-DRS SQ-SD 
中 华 蜜蜂 刚 出 房 .性 成 熟 蜂王 和 雄 蜂 之 间 差 异 表达 基因 的 数目 
Fig.3 Number of differentially expressed genes among newly emerged and sexually matured drones and queens of Apis cerana cerana 
ED: 刚 出 房 雄 蜂 Newly emerged drone; EQ: 刚 出 房 蜂王 Newly emerged queen; SD: 性 成 熟 雄 蜂 Sexually matured drone; DRS: 雄 蜂 生殖 系统 


SD-DRS 








Drone reproductive system; QRS: 蜂王 生殖 系统 Queen reproductive system; SQ: 性 成 熟 蜂 王 Sexually matured queen. 


比较 性 成 熟 蜂 王 与 刚 出 房 蜂王 发 现 , 两 时 期 存 
在 14 072 个 差异 基因 ,其 中 11 188 个 基因 在 性 成 熟 
蜂王 中 上 调 ,2 884 个 基因 下 调 , 包 括 27 个 表皮 有 蛋 
白 7 个 热 应 激 和 蛋白 ,6 个 气味 结合 蛋白 和 至 少 15 个 
细胞 色素 P450 等 基因 。 人 性 成 熟 雄 蜂 与 刚 出 房 雄 蜂 
相 比较 ,有 9 741 个 基因 存在 差异 表达 ,其 中 6 405 
个 基因 在 性 成 熟 期 上 调 表 达 ,3 336 个 基因 下 调 表 
达 , 包 括 14 个 表皮 蛋白 .9 个 气味 结合 蛋白 和 至 少 9 
个 热 应 激 蛋 白 、10 个 细胞 色素 P450。 

刚 出 房 的 蜂王 与 刚 出 房 的 雄 蜂 相 比 ,它们 之 间 有 
8 331 个 基因 存在 差异 表达 ,其 中 在 刚 出 房 蜂王 中 上 
调 表达 的 基因 有 2 406 个 ,在 刚 出 房 雄 蜂 中 上 调 表达 
的 基因 有 5 925 个 。 然 而 ,到 性 成 熟 期 之 后 ,两 者 存在 
15 487 个 差异 基因 ,在 蜂王 和 雄 蜂 中 上 调 表达 的 基因 
数量 分 别 为 7 865 和 7 622 个 。 这 些 差异 基因 中 包括 
了 9 种 气味 结合 蛋白 ,其 中 只 有 气味 结合 蛋白 基因 15 
(0BP15 ) 在 性 成 熟 蜂王 中 上 调 表 达 , 其 余 8 种 气味 结 
合 蛋 白 相 关 基 因 (0BP9, OBP1, 0BP6, 0BP14， 
0BP21 ,ASP6,ASP4 和 4SP5 ) 在 性 成 熟 雄 蜂 中 上 调 
表达 。 

研究 发 现 ,性 成 熟 期 蜂王 与 雄 蜂 生殖 系统 之 间 
存在 26 976 个 差异 表达 基因 ,其 中 ,在 蜂王 生殖 系 
统 中 高 表达 的 基因 有 14 554 个 ,另外 12 422 个 基因 
在 雄 蜂 生殖 系统 中 高 表达 ,其 比例 与 性 成 熟 期 蜂王 
和 雄 蜂 差异 基因 比例 相似 。 在 这 些 差异 基因 中 ,有 









































7 个 气味 结合 蛋白 存在 显著 差异 表达 ,其 中 只 有 
0BP14 基因 在 蜂王 生殖 系统 中 上 调 表 达 , 另 外 6 个 
气味 结合 相关 基因 (OBP9,OBP1,OBP16, 0BP21， 
ASP2 和 ASP4 ) 均 在 雄 蜂 生殖 系统 中 上 调 表 达 。 

2.4 GO 和 KEGG 富 集 分 析 

将 差异 基因 作为 前 景 ,全 部 基因 作为 背景 ,进行 
GO 和 KEGG 的 富 集 分 析 。 性 成 熟 蜂 王 与 刚 出 房 蜂 
王 差 异 基 因 按 照 细 胞 组 成 .分 子 功能 和 生化 过 程 三 
大 类 , 共 分 为 56 个 功能 群 ,分 别 有 28,10 和 50 个 
GO term 在 三 大 类 中 显著 富 集 (P <0.05) ,包括 气味 
结合 类 。 性 成 熟 雄 蜂 与 刚 出 房 雄 蜂 差异 基因 按照 生 
物 合成 .细胞 组 成 和 分 子 功能 三 大 类 , 共 分 为 55 个 
功能 群 , 有 26,23 和 50 个 GO term 在 三 大 类 中 显著 
富 集 (P<0.05), 包 括 气 味 结合 类 。 人 性 成 熟 蜂 王 与 
性 成 熟 雄 蜂 差异 基因 按 生物 合成 .细胞 组 成 和 分 子 
功能 三 大 类 分 为 56 个 功能 群 ,有 13,6 和 29 个 CO 
term 在 三 大 类 中 显著 富 集 (已 <0.05) ,包括 气味 结 
合 类 。 性 成 熟 期 的 蜂王 生殖 系统 与 雄 蜂 生 殖 系统 差 
异 基因 共 分 57 个 功能 群 ,有 18,15 和 46 个 GO term 
在 三 大 类 中 显著 富 集 (P <0.05)。 

将 差异 基因 与 KEGG 数据 库 进 行 比 对 ,性 成 熟 
蜂王 与 刚 出 房 蜂王 之 间 的 差异 基因 能 归 类 到 255 个 
生化 途径 ,其 中 40 个 显著 富 集 (P <0.05); 性 成 熟 
雄 蜂 与 刚 出 房 雄 蜂 之 间 的 差异 基因 能 归 类 到 257 个 
生化 途径 ,其 中 67 个 显著 富 集 (P < 0.05) ;性 成 熟 
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蜂王 与 性 成 熟 雄 蜂 差 异 基因 能 归 类 到 256 个 生化 途 
径 , 其 中 51 个 显著 富 集 (P <0.05); 性 成 熟 期 的 蜂 
王 生 殖 系 统 与 性 成 熟 期 的 雄 蜂 生殖 系统 差异 基因 能 
归 类 到 257 个 生化 途径 ,其 中 49 个 显著 宣 集 (P < 
0.05) 。 
2.5 ”qRT-PCR 验证 结果 

从 差异 基因 中 选择 了 4 种 热 应 激 和 蛋白 基因 和 
TPX 基因 在 性 成 熟 蜂王 与 刚 出 房 蜂王 以 及 性 成 熟 蜂 





























王 生殖 系统 与 性 成 熟 雄 峰 生 殖 系统 样品 中 进行 验 
证 。4 个 热 应 激 和 蛋白 基因 和 TPX 基因 在 性 成 熟 蜂王 
中 的 表达 均 显 著 高 于 刚 出 房 蜂王 ( 表 5) ,与 转录 组 
结果 一 致 。 另 外 ,HSP70Ab 和 HSP23 基因 在 性 成 熟 
雄 蜂 生殖 系统 显著 上 调 表 达 , 而 HSP90 和 TPX 基因 
在 性 成 熟 蜂王 生殖 系统 中 显著 上 调 表达 ( 表 6) ,与 
转录 组 结果 一 致 。HSP83 基因 在 性 成 熟 的 两 种 生殖 
系统 表达 量 差异 不 显著 ( 表 6) 。 

















表 5 5 个 基因 在 中 华 蜜蜂 蜂王 中 的 表达 


TableS Quantitative analysis of five genes in Apis cerana cerana queens 


































































































刚 出 房 蜂王 转录 组 RPKM 值 ”性 成 熟 蜂王 转录 组 RPKM 值 ” 刚 出 房 蜂 王 基因 表达 量 。 性 成 熟 蜂王 基因 表达 量 , 
基因 2 1 2 P 值 
RPKM of transcriptome in RPKM of transcriptome in Expression level of gene in Expression level of gene in 
Gene P-value 
newly emerged queen sexually matured queen newly emerged queen sexually matured queen 
HSP70Ab 5.6424 20. 6883 1 3.6106 <0.05 
HSP90 24.7959 78.4251 1 3.2917 <0.05 
HSP83 70.4126 304.5616 1 9.0204 <0.05 
TPX 941.1732 2 113.4686 1 4.6898 <0.05 
HSP23 123.0106 922.3816 1 4.1281 <0.05 
表 6 5 个 基因 在 中 华 蜜蜂 性 成 熟 蜂 王 和 雄 蜂 生殖 系统 中 的 表达 
Table6 Quantitative analysis of five genes in the sexually matured reproductive system 
of Apis cerana cerana queens and drones 
其 大 雄 蜂 生殖 系统 转录 组 RPKM 值 蜂王 生殖 系统 转录 组 RPKM 值 雄 蜂 生殖 系统 基因 表达 量 蜂王 生殖 系统 基因 表达 量  P 值 
C RPKM of transcriptome in RPKM of transcriptome in Expression level of gene in Expression level of gene in Pival 
ek drone reproductive system queen reproductive system drone reproductive system queen reproductive system A 
HSP70Ab 16. 8763 0. 4868 1 0. 0338 <0.05 
HSP90 164. 6236 501. 891 1 1.4132 <0.05 
HSP83 83. 5527 406. 2492 1 1. 6983 >0.05 
TPX 971. 6687 3 498. 111 1 1.7082 <0.05 
HSP23 871.4282 40. 1814 1 0.0158 <0.05 
素 具 有 必然 关系 ,从 而 有 利于 蜜蜂 正常 交尾 。 
全， My 
3 讨论 昆虫 表皮 蛋白 是 结构 蛋白 , 它 和 几 丁 质 一 起 组 








本 研究 通过 高 通 量 测序 技术 对 中 华 蜜蜂 蜂王 和 
雄 蜂 性 成 熟 时 期 及 性 成 熟 期 生殖 系统 的 转录 组 差异 
进行 了 分 析 。 测 序 结 果 发 现 :所 有 试验 样本 高 质量 
的 数据 ( Q20) 均 高 于 94% ,说 明 测 序 质量 较 好 。 所 
有 转录 组 测序 获得 的 reads ,组 装 后 ,获取 90 839 个 
unigenes ,平均 长 度 1 549 bp, 而 且 大 部 分 unigenes 
能 比 对 到 亲缘 关系 较 近 的 物种 , 仅 7.6% unigenes 
未 知 ,这 也 说 明 测序 组 装 质 量 较 好 。 

对 差异 表达 基因 进行 分 析 发 现 ,性 成 熟 的 蜂王 
和 雄 蜂 与 刚 出 房 时 期 相 比 , 均 存 在 大 量 显著 表 达 差 
异 的 基因 ,而 且 表 皮蛋 白 、 热 应 激 和 蛋白、 气味 结合 蛋 
白 以 及 细胞 色素 P450 基因 均 存 在 差异 表达 ,并 且 气 
味 结合 类 显著 富 集 ,说 明 这 些 基因 与 雄 蜂 识别 信息 







































































成 昆虫 抵御 外 界 环 境 的 屏障 ,形成 角质 层 ( 刘 清明 
等 ,2010)。 本 研究 发 现在 雄 蜂 或 蜂王 生长 发 育 过 
程 中 差异 表达 的 表皮 蛋白 基因 均 在 刚 出 房 时 期 上 调 
表达 ,主要 原因 在 于 刚 出 房 的 蜜蜂 比较 柔弱 ,体格 较 
软 , 机 体 有 待 于 进一步 发 育成 熟 ,从 而 需要 表皮 蛋白 
基因 表达 活跃 ,基因 表达 产生 的 蛋白 经 过 转运 机 制 
运输 到 昆虫 体 辟 上皮 细胞 之 外 而 参与 外 骨骼 的 
构建 。 

细胞 色素 P450s 广泛 分 布 于 各 种 需 氧 生物 中 ， 
如 昆虫 植物 细菌 等 。 它 在 昆虫 中 除了 参与 体内 的 
生物 代谢 外 ,还 具有 激素 的 合成 与 分 解 ( 李 添 等 ， 
2004) 。 细 胞 色素 P450 基因 大 部 分 在 蜜蜂 性 成 熟 
期 高 度 表 达 , 可 能 参与 了 蜜蜂 体 表 信 息 素 的 合成 ,如 
9- 氧 代 - 癸 二 烯 酸 (9-0ODA)。 男 外 , 细胞 色素 
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CYP342A1,CYP6BD1 和 CYP6AS5 基因 在 性 成 熟 期 
的 蜂王 和 雄 蜂 均 上 调 表 达 , 而 细胞 色素 CYP4AV1 
基因 却 下 调 表 达 , 这 些 基 因 可 能 是 蜜蜂 性 成 熟 期 的 
部 分 管家 基因 。 细 胞 色素 P450 还 具有 抗 氧 化 和 加 
强 免疫 的 作用 ,降低 外 来 物 的 侵害 ,还 可 能 与 热 应 激 
和 蛋白 一 起 抵抗 外 界 不 良 环境 的 影响 。 

气味 结合 蛋白 是 一 类 识别 或 结合 外 界 气 味 分 子 
并 传输 给 受 体 的 一 类 蛋白 ,广泛 分 布 于 昆虫 触角 器 
官 或 表皮 细胞 中 。 本 研究 发 现 ,大 部 分 气味 结合 蛋 
白 基 因 在 性 成 熟 雄 蜂 上 调 表 达 , 这 可 能 与 蜜蜂 交尾 
机 制 有 关 , 雄 蜂 为 了 在 空中 有 效 地 追逐 蜂王 并 进行 
交尾 ,就 需要 获取 蜂王 在 婚 飞 过 程 中 释放 的 信息 素 。 
有 些 上 调 基因 或 下 调 基因 是 其 不 同时 期 的 特异 性 表 
达 ( 李 红 亮 ,2007 ) 。 

性 成 熟 时 期 是 蜂王 和 雄 蜂 成 功 交 配 的 必然 条 件 
之 一 ,而 且 其 体 表 会 释放 信息 素来 吸引 对 方 。 比 较 
性 成 熟 期 的 蜂王 和 雄 蜂 差异 表达 基因 发 现 ,两 者 之 
间 存 在 15 487 个 差异 基因 ,包括 9 种 气味 结合 蛋白 
基因 ,其 中 8 种 在 雄 蜂 中 上 调 表 达 , 说 明雄 蜂 性 成 熟 
期 嗅觉 特别 敏感 ,有 利于 婚 飞 时 进行 竞争 交配 。 同 
样 ,比较 两 者 之 间 的 生殖 系统 发 现 , 有 7 种 气味 结合 
蛋白 基因 存在 差异 表达 ,其 中 6 种 在 雄 蜂 生殖 系统 
上 调 表达 ,这 也 说 明 气 味 结合 蛋白 与 雄 蜂 性 腺 发 育 
有 着 密切 的 关联 。 

另外 , 刚 出 房 雄 蜂 与 刚 出 房 蜂 王 存在 9 331 个 
差异 表达 基因 ,而且 其 中 5 925 个 基因 在 刚 出 房 雄 
蜂 中 上 调 表达 ,这 可 能 是 由 于 刚 出 房 雄 蜂 新 陈 代 谢 
比 刚 出 房 蜂王 要 快 ,这 有 待 于 进一步 研究 。 

本 研究 利用 高 通 量 测序 技术 对 中 华 蜜蜂 蜂王 和 
雄 蜂 在 刚 出 房 .性 成 熟 时 期 以 及 其 性 成 熟 期 时 期 生 
殖 系 统 的 转录 组 进行 测序 和 分 析 , 揭 示 了 中 华 蜜蜂 
发 育 的 两 个 关键 时 期 和 性 成 熟 时 期 生殖 系统 转录 组 
的 整体 表达 特征 。 初 步 发 现 中 华 蜜蜂 在 发 育 过 程 
中 ,其 表皮 蛋白 基因 、 热 应 激 蛋 白 基 因 、 细 胞 色素 
P450 基因 以 及 气味 结合 和 蛋白 基因 表达 发 生 了 明显 
的 变化 ,为 深入 研究 中 华 蜜蜂 生长 发 育 、 繁 殖 以 及 相 
关 基 因 功 能 及 调控 等 提供 了 丰富 的 数据 资源 。 
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